Nonlocal interactions in doped cuprates: correlated motion of Zhang— Rice polarons 
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In-plane, inter-carrier correlations in hole doped cuprates are investigated by ab initio multiconfig- 
uration calculations. The dressed carriers display features that are reminiscent of both Zhang-Rice 
(ZR) CuC>4 singlet states and Jahn-Teller polarons. The interaction between these quasiparticles 
is repulsive. At doping levels that are high enough, the interplay between long-range unscreened 
Coulomb interactions and long-range phase coherence among the O-ion half-breathing vibrations on 
the ZR plaquettes may lead to a strong reduction of the effective adiabatic energy barrier associated 
to each polaronic state. Tunneling effects cannot be neglected for a relatively flat, multi-well energy 
landscape. We suggest that the coherent, superconducting quantum state is the result of such co- 
herent quantum lattice fluctuations involving the in-plane O ions. Our findings appear to support 
models where the superconductivity is related to a lowering of the in-plane kinetic energy. 



The parent compounds of the high-temperature 
cuprate superconductors are Mott insulating antiferro- 
magnets. Superconductivity (SC) occurs in these sys- 
tems upon hole or electron doping. While the mag- 
netic interactions are remarkably strong, there is increas- 
ing evidence that many properties cannot be understood 
without taking into account the lattice degrees of free- 
dom. Experiments that indicate strong elect ron-phonon 
(EP) couplings are the inelastic neutron scattering 0,0, 
angle-resolved photoemission (PE) 0, Q , X-ray absorp- 
tion fine structure (XAFS) J5j, and electron paramag- 
netic resonance |6:j. Anomalies in the phonon, PE, and 
XAFS spectra were previously addressed with t—J or 
Hubbard-like models supplemented with EP interaction 
terms 0,11,11,^3. Strong and anomalous electron-lattice 
couplings were also found by ab initio, explicitly corre- 
lated calculations [TlllT^ |. In the case of hole doping [lflj . 
multiconfiguration (MC) calculations on clusters of few 
CuC>6 octahedra show that the doped holes enter O 2p x 
and 2p y orbitals that form cr-bonds with the open-shell 
Cu id x 2_ y 2 orbitals [ill ] and give rise to singlet states 
similar to the kind of confi gura tion proposed by Zhang 
and Rice (ZR) long time ago 14]. A major deviation from 
the original ZR picture is that the formation of such a 
singlet is associated with significant lattice deformations. 
The most stable configuration corresponds to Cu-0 dis- 
tances that are shorter by 5-6%. Most remarkably, O-ion 
displacements that restore the translational symmetry in- 
duce strong charge redistribution. The 2p hole, which for 
a distorted Cu04 plaquette is equally distributed over 
the four anions, can be partially transferred onto a sin- 
gle ligand to give an electronic wavefunction (WF) with 
a dominant contribution from a ...-Cu<i 9 -Op 5 -Cu<i 9 - 
... configuration. The ZR-like singlet polaron (ZRP) can 
hop thus within the Cu02 plane via such cP-p 5 - d 9 states 
through coupling to the oxygen vibrations. The energy 
barrier associated with the hopping process is few hun- 
dreds meV, for an isolated 2p hole \Vi\ . 

d 9 -p 5 -d 9 entities were postulated to exist in cuprates 
by several authors [TsL Hq . This single-oxygen 2p hole 
and the adjacent Cu holes are viewed as a spin-1/2 quasi- 



particle state in which the 3d spins are ferromagnetically 
(FM) polarized, i.e. coupled to a triplet, and the spin 
on the O ion is low-spin coupled to it. Interactions that 
could lead to pairing of these quasiparticles and SC were 
analyzed in 0, 0] in terms of p -d, extended Hubbard 
Hamiltonians. It was pointed out [iall3 that such mod- 
els imply different low-energy physics as compared to 
the effective single-band model proposed by ZR [lij]. In- 
terestingly, the ab initio calculations 01 snow that the 
two concepts, the ZR and the d 9 -p 5 -d 9 configurations, 
are in fact intimately related. The connecting element 
is the lattice degree of freedom, more precisely the Cu- 
O bond-length fluctuations. Near-degeneracy effects be- 
tween a quadrisinglet (QS) polaronic state similar to the 
ZR state and a two-site spin-singlet polaron were stud- 
ied in the frame of a Holstein-Hubbard model in [l8| . 
Although some details are quite different, the two-site 
spin-singlet resembles a two-hole d 9 -p 5 singlet configu- 
ration. It was shown that in the parameter region with 
near-degeneracy, the effective mass of the QS polaron is 
sharply reduced. Near-degeneracy plus quantum tunnel- 
ing effects associated with a rather fiat energy landscape 
are believed to be 0] crucial ingredients in the theory 
of "bipolaronic" SC. 

In this Letter we investigate interactions between two 
or more spatially separated ZRPs. We find that long- 
range phase coherence among the O-ion half-breathing 
vibrations on the ZR CUO4 plaquettes may lead to a 
correlated motion of the 2p holes. Several effects come 
here into play. The local electron-lattice couplings make 
possible the transfer of the ZR singlet to an adjacent pla- 
quette. The longer-range inter-carrier interactions are re- 
sponsible, at dopant concentrations that are high enough 
and in a "dynamic" regime, for a significant lowering of 
the effective energy barrier seen by each ZR quasiparticle. 
This is exactly the kind of picture proposed, in a some- 
what different context, by Hirsch [19|: single carriers are 
heavily dressed at very low concentrations by the interac- 
tion with the local environment; they partially "undress" 
and their effective mass decreases at higher concentra- 
tions due to collective inter-carrier interactions. 
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Our analysis is mainly based on results obtained by 
ab initio electronic structure methods from traditional 
quantum chemistry. The calculations are performed on 
finite clusters. The many-electron WF is constructed as a 
full configuration-interaction (CI) expansion in the space 
defined by a limited set of so-called active orbitals. These 
are chosen as those orbitals that are expected to con- 
tribute to degeneracy effects, that is, strong mixing be- 
tween configurations which have the same, or nearly the 
same, energy The rest of the orbitals are either doubly 
occupied in all configurations or empty. The former set of 
orbitals is called inactive, the latter is referred to as the 
virtual orbital set and spans the rest of the orbital space, 
defined from the basis set used to build the "molecular" 
orbitals (MOs). This is known as the Complete Active 
Space (CAS) MC approach [2(j . The variational param- 
eters of the WF are both the CI and MO coefficients. 

Inter-carrier interactions are studied by doping a 7- 
plaquette [CU7O22] linear cluster with two holes. The 
length of the cluster corresponds roughly to the scale of 
the superconducting coherence length. We used the lat- 
tice parameters of the hole doped Lai.85Sro.isCu04 com- 
pound |2l|. with the in-plane lattice constant a = 3.78 A. 
The 7-plaquette cluster is embedded in an array of for- 
mal ionic point charges (PCs) at the experimental lat- 
tice positions that reproduce the crystal Madelung po- 
tential. The Cu 2+ and La 3+ nearest neighbors (n.n.) 
are represented by effective model potentials |22f. Trial 
calculations were first performed with embedding PCs 
that correspond to a fully ionic picture of the undoped 
material, La2Cu04. In this situation, due to the large 
mutual repulsion, the CAS Self-Consistent Field (SCF) 
calculations converge to solutions where the 2p holes are 
each located on the outermost CUO4 plaquettes. How- 
ever, at not too small doping other 2p holes exist in the 
immediate vicinity. In a next step, we decided to take 
into account the presence of other doped holes by intro- 
ducing few extra elementary positive charges around the 
cluster. The results discussed in the following paragraphs 
were obtained with a PC embedding where eight positive 
charges were added. The positions of these extra PCs are 
shown in Fig. 1. The reason behind choosing this type of 
arrangement is discussed below. 

It is widely accepted j2^| that charge and magnetic 
inhomogeneities appear in cuprates at least dynamically. 
Goodenough [25| argued that at optimal doping (x~ 1/6) 
and low temperatures, 1/3-doped (fluctuating) stripes 
are formed on every second chain of plaquettes. In a 
ZR-like picture, the 1/3-doped stripes could be repre- 
sented as in Fig. 1(a), for example. The SC is associated 
in Ref. 251 with the formation of traveling charge-density- 
waves [24| along the doped chains. In models with EP 
coupling, the relevant phonon modes would be the O-ion 
half- breathing modes with q = (1/3 , , 0). Experimental 
data that appear to support Goodenough's model are 
the strong anomalies at such wave vectors, q=(0.25- 
0.33,0,0) 0,0- These anomalies are actually observed 
in both hole and electron doped systems Q- Coming 




FIG. 1: (a) 1/3-doped stripes at optimal doping in a ZR- 
like picture, see text, (b) The 7-plaquette cluster employed 
for the CAS SCF calculations and the positions of the extra 
positive charges (the filled circles) within the PC embedding. 
These correspond to the situation where every second pair of 
ZRPs along the 1/3-doped stripes is removed. One of the two 
symmetry-equivalent lowest-energy configurations is sketched 
for the "quantum" cluster. Only the Cu sites are displayed. 
The squares represent distorted ZR plaquettes. 



back to the positive charges added to our embedding, 
they would correspond to the situation where every sec- 
ond pair of ZRPs in Fig. 1(a) is removed, i.e. x = 1/12 
doping. In addition, we retain only the n.n. pairs around 
the "quantum" cluster. For symmetry reasons, the extra 
positive charges on the left-hand side are shifted by one 
plaquette to the right. 

With the extra positive charge around the cluster, the 
two O holes are pushed from the outermost plaquettes 
toward the inner region. To study electron-lattice cou- 
plings, we first determined the cluster geometry config- 
uration that minimizes the total energy. However, only 
the positions of those ions are relaxed that have each 
neighbor represented at the all-electron level. These are 
the Cu ions and the intervening ligands O&. Under the 
assumption of charge segregation into hole-rich and hole- 
poor regions and filamentary conduction [Tol llol l23l |25| 
the vibrations of these ions should be the most relevant. 
The CAS SCF calculations were carried out with a mini- 
mal active space. That is one orbital per hole, in our case 
nine electrons in nine orbitals. It was shown 12] that the 
overall picture does not change when using larger active 
spaces. Extra orbitals added to the active space turn out 
to have occupation numbers either very close to 2.00 or 
to 01 • Also, when enlarging the active space, it is 
not obvious how many extra orbitals are required for a 
balanced description of the different configurations. All 
calculations were performed with the MOLCAS 6 package 
|26|. All-electron atomic natural-orbital basis sets from 
the MOLCAS library were employed, with the following 
contractions: Cu 21sl5pl0d/5sAp3d, O 14s9p/4s3p. As 
already mentioned, the apical oxygens are not included 
in the "quantum" cluster, but modeled by formal PCs. 
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TABLE I: MPs illustrating the distribution of the O-holes. 
The lowest-energy configuration corresponds to distorted, ZR 
plaquettes separated by two other Cu 3d*2_ y a ions (1st col- 
umn). The ZR plaquettes are labeled with indices 1 and 2. 
Other notations are as in Fig. 1. When the Oj, ions are shifted 
back to the middle positions, each hole is partially transferred 
onto a single ligand (2nd column). In parentheses, the same 
quantities are shown for a single doped hole, see text. 



Configuration 


Dist. 


Cu-O b , 


Undist. str. 




ZRP states 




Cui 3d x 2_ y 2 


1.03 


(1.03) 


1.09 (1.08) 


0[ b 2 Py 


1.64 


(1.62) 


1.77 (1.76) 


Ola 2p x (x2) 


1.61 


(1.62) 


1.69 (1.68) 


OJi 2 Py 


1.61 


(1.62) 


1.29 (1.31) 


CU2 3d x 2_y2 


1.03 




1.09 


O l 2b 2 Py 


1.65 




1.77 


2a 2p x (x2) 


1.61 




1.68 


05 & 2 P y 


1.60 




1.30 


Rel. En. (meV) 


0.00 


(0.00) 


650 (380) 



The CAS SCF calculations always converge to states 
where the 2p holes are separated by more than two lat- 
tice constants. However, depending on geometry and the 
initial guess for the orbitals, these states may have dif- 
ferent character. Two symmetry-equivalent minimum- 
energy geometries were identified. These are configura- 
tions where two ZRP states are formed on Cu04 pla- 
quettes that are separated by two other Cu ions, see 
Fig. 1(b). The spins on these intervening cations are 
antifcrromagnetically (AFM) coupled. For the ZR pla- 
quettes, the Cu-O b distances are shorter by 5%. Mul- 
liken populations (MPs) of the relevant Cu 3d and O 2p 
atomic orbitals on these plaquettes are listed in Table I 
(first column). The composition of the bonding (B) and 
antibonding (AB) d x 2_ y -2—p xy er-orbitals is depicted in 
Fig. 2. Due to correlation effects, the occupation num- 
bers (ONs) of the AB orbitals are relatively large. Since 
the charge distribution on the two ZR plaquettes is nearly 
identical, only the orbitals on one of these plaquettes are 
shown in the figure. 

When the bridging Ob ions arc shifted back to the 
middle positions [Tj, each of the 2p holes is partially 
transferred onto a single ligand. This is illustrated in 
Table I and Fig. 2. The results were obtained by static 
total energy calculations (adiabatic approximation), ft 
can be seen in Fig. 2 that the doped-hole orbitals have 
now much weight onto one of the adjacent plaquettes as 
well. O-ion half-breathing displacements that shorten the 
Cu-Oj, bonds on these adjacent plaquettes induce prop- 
agation of the ZR singlet (s) along the Cu-0&-... chain. 
We found that the ZR 2p hole is already moved to the 
other plaquette for distortions of less than 0.5% from the 
high-symmetry geometry. The adiabatic energy barrier 
is in the two-hole cluster 325 meV per hole, see Table I. 
Separate calculations were performed with a single doped 
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FIG. 2: (a)-(b) B and AB orbitals defining a ZR-like state on 
a distorted plaquette. Their ONs are 1.86 and 0.14. (c)-(d) 
B and AB orbitals in the undistorted structure. The ONs are 
1.74 and 0.26. The 2p hole has now the largest weight onto a 
single "chain" oxygen, see also Table I. The whole 7-plaquette 
cluster is drawn. The phases of the atomic orbitals are shown 
with different colors. 



hole and no extra positive charge in the PC embedding. 
The energy barrier is in this case 380 meV [2jj . The en- 
ergy difference between the two situations, 55meV per 
hole, is considerable. In the presence of other mobile 
2p hole polarons in the Cu02 plane(s), the effect would 
be further enhanced and should lead to a strong reduc- 
tion of the polaron effective mass. We arrive thus at 
Hirsch's paradigm 0, |2^|: carriers are not happy be- 
ing together, but they are happy moving together. This 
was described as a kinetic-energy-driven mechanism for 
SC and a change of color was predicted when entering 
the superconducting state [l^- That the onset of SC 
is indeed accompanied by a transfer of spectral weight 
from the visible region to the infrared range was recently 
confirmed by Molegraaf et al. |2£j . 

In comparison to the study from Ref.|l9|| and other 
model Hamiltonian studies, we are able to identify much 
more precisely the crucial interactions in the doped plane. 
The CAS SCF results give also us hints about the nature 
of the pairing mechanism. It is obvious that this is re- 
lated to both long-range Coulomb repulsion effects and 
shorter-range interactions. Our results may be associated 
with the 4-plaquette correlation bag of Goodenough [2£j . 
The ZR singlet polarons in such a correlation bag are 
separated by AFM coupled d x 2_ y 2 spins that bear some 
resemblance to the spin-singlet pairs of the resonating- 
valence-bond state 30]. Cu-0 bond- length fluctuations 
may induce charge transfer along a chain of plaquettes. 
In the intermediate, undistorted geometry the quasilo- 
calized 2p holes tend to polarize FM the n.n. 3d spins. 
In other words, the 3c? spin on the adjacent plaquette is 
already "prepared" for the transfer of the 2p hole. Also, 
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strong antiferromagnetic correlations exist between the 
d 9 -p 5 -d 9 units and the "bipolaronic" WF has a complex 
structure. Our data suggest that two ZR-like quasipar- 
ticles could form a 4-plaquette bipolaronic singlet state 
only above a certain hole concentration. The distribu- 
tion of "doped" holes displayed in Fig. 1(b) corresponds 
(at the scale of the figure) to 8-9% doping, only slightly 
higher than the concentration(s) where SC actually oc- 
curs. We also note that the effects described above, i.e. 
a lowering of the effective energy barrier due to longer- 
range interactions, are still obtained when the positions 
of the extra positive embedding PCs are changed, indi- 
vidually or collectively, by distances as large as one lattice 
constant. 

More investigation is needed for a better understand- 
ing of these systems and quantitative estimates. First, a 
more rigorous treatment of the dynamical electron cor- 
relation [23| is required. Secondly, the electronic and 
lattice degrees of freedom should be described on equal 
footing, beyond the adiabatic approximation. Another 
non-trivial task is a self-consistent procedure to deter- 
mine the (dynamic) in-plane distribution of the charge 
carriers, pe rhap s in the spirit of the dynamical mean- 
field theory |31| . Nevertheless, the purpose of the present 
study is not to provide highly accurate numbers, but to 
gain better insight into the very nature of the many-body 
in-plane interactions. It is our belief that realistic mod- 
els can only be constructed with the help of this type of 
first- principles investigations. 



To conclude, using a WF-based ab initio approach, we 
are able to describe microscopic electron-lattice, mag- 
netic, and nonlocal Coulomb interactions characteriz- 
ing the hole dynamics in (under)doped cuprates. The 
dressed carriers display features which are reminiscent 
of both ZR spin-singlet states and Jahn-Teller polarons. 
We find that at doping levels that are large enough, 
the adiabatic energy barrier associated to each polaronic 
state is significantly lower for the synchronized, coherent 
"hopping" of two ZRPs than in the case of an isolated 
ZRP. This is due to unscreened inter-carrier Coulomb 
interactions and requires phase coherence among the O- 
ion half-breathing vibrations. Tunneling effects cannot 
be neglected for a relatively flat, multi-well energy land- 
scape 13 • Following the discussion from Ref . [Tsj] . 
we suggest that the coherent, superconducting phase is 
the result of such coherent quantum lattice fluctuations 
involving the in-plane O ions. Coherent, collective tun- 
neling effects in the Cu02 planes should lead to extended, 
"resonant" states and a lowe ring of the in-plane ki- 
netic energy, see also 0, |2^, EIj ■ Polaronic behav- 
ior and local double-well potentials were also found in 
the electron doped case 1 1 *Ji . which suggests some com- 
mon "global" features for the two types of doping: quan- 
tum polarons coupled through spin and longer-range un- 
screened Coulomb interactions. 

We thank I. Eremin and P. Fulde for fruitful discus- 
sions. 
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